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The mechanism of rapid and selective heating of magnetic metal oxides under the magnetic 
field of microwaves which continues beyond the Curie temperature T c is identified by using the 
Heisenberg model. Monte Carlo calculations based on the energy principle show that such heating 
is caused by non-resonant response of electron spins in the unfilled 3d shell to the wave magnetic 
field. Small spin reorientation thus generated leads to a large internal energy change through the 
exchange interactions between spins, which becomes maximal around T c for magnetite Fe30<i. The 
dissipative spin dynamics simulation yields the imaginary part of the magnetic susceptibility, which 
becomes largest around T c and for the microwave frequency around 2 GHz. Hematite Fe203 with 
weak spontaneous magnetization responds much less to microwaves as observed in experiments. The 
heating of titanium oxide by microwave magnetic field only when oxygen defects are present is also 
explained by our theory in terms of the absence of spontaneous magnetization. 
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I. INTRODUCTION 

Microwave sintering is the process in which electro- 
magnetic energy of microwaves is delivered directly to 
electrons or atoms in materials. It enables rapid heat- 
ing that leads to melting and recrystallization of original 
solid matters. A laboratory experiment demonstrated 
that a grained magnetite sample was heated to above 
1300°C within a minute by applied microwaves^. The 
amount of energy consumption can be reduced and the 
release of CO2 gas be halved in the microwave iron mak- 
ing which reduces energy and environmental problems in 
modern iron industry^. Later it was shown that various 
metallic oxides including magnetite and titanium oxide 
with oxygen defects Ti02- X {x > 0) were sintered quickly 
at the magnetic field maximum (i.e. the electric field 
null) in the microwave cavity experiments 3 -. The prod- 
uct sintered in the microwave magnetic field was shown 
to have good magnetic properties as magnets. 

The sintering of magnetite/hematite powders of mi- 
cron sizes is characterized by selective heating. Only the 
domain of magnetite FeaCU with strong magnetization 
was heated to 1300°C in the microwave magnetic field 4 , 
which is much above the Curie temperature T c (585°C). 
The adjacent domain of hematite Fe2C>3 remained at low 
temperatures. The microwave sintering of magnetic ma- 
terials is a non-resonant process that occurs with or with- 
out a static magnetic field and for any amplitude of wave 
magnetic field at a fixed microwave frequency. With fer- 
romagnetic metals, microwave resonance absorption^ and 
the non-resonant bolometric effect on dc resistances^ 
were observed in ferromagnetic resonance (FMR) exper- 
iments under microwaves. These phenomena required a 
static magnetic field of specific strength to have reso- 
nance with microwaves, and were observed mainly below 
room temperature, which are attributed to Joule heating 
or eddy current. The microwave sintering, which occurs 



under wide range of strength and frequency of microwave 
magnetic field without a static magnetic field is physi- 
cally a different process from the FMR process. 

To date, the heating of dielectric materials by the mi- 
crowave electric field was extensively studied^. We ex- 
amined the heating of water, salt-contained water and ice 
by microwaves using molecular dynamics simulatio n 10 ' 11 . 
The heating was attributed, respectively, to excitation 
of electric dipoles, the Joule heating of salt ions, and 
the weakening of the hydrogen-bonded H2O network by 
the presence of salt ions. (Pure ice was not heated by 
2.5GHz microwave due to tight hydrogen bonds of wa- 
ter molecules.) For bulk metals, microwaves are reflected 
at the surface due to skin effects and do not transfer 
energy, while they penetrate into grained metallic pow- 
ders (10/^m in diameter) for a few centimeters and heat 
themi 2 -. However, the sintering mechanism of metal 
oxides that have spontaneous magnetization under the 
magnetic field of microwaves has not been resolved. 

In this paper, we show theoretically that the sintering 
of particles of magnetic metal oxides, including magnetite 
and titanium oxide with oxygen defects, by a microwave 
magnetic field is due to non-resonant response of electron 
spins in the unfilled 3d shell. We first use the energy prin- 
ciple and perform the Monte Carlo simulations. A small 
spin perturbation in response to an alternating external 
magnetic field results in a large internal energy change 
through exchange interactions. Next, studies in the time 
domain are done by performing dissipative spin dynam- 
ics simulation and detecting the linear response of spins. 
The temperature and frequency dependence of the heat- 
ing rate by the microwave magnetic field is obtained on 
the basis of the imaginary part of magnetic susceptibil- 
ity Xi- These results agree well with those by the energy 
principle, and Xi quantitatively accounts for the rapid 
heating of magnetite in the sintering experiments. 



II. NUMERICAL PROCEDURES 



We have used the following procedures in our nu- 
merical simulations. The magnetization of magnetite 
and hematite is well described by the Heisenberg model 
above the Verwey transition temperature (120 K)^2 since 
electrons are roughly localized^. The internal energy 
U of the magnetic system is represented by the three- 
dimensional spin vector of the electron at the i-th site, 
the exchange interaction coefficient Jy between the i-th 
and j-th sites, and the external magnetic field B m of mi- 
crowaves, which reads 



U (B w ) = — y] JijSi ■s j + } j guBSi ■ B t 



(1) 



The summation of the exchange interactions in the first 
term is taken over the pairs of nearest neighbor sites. 
The magnitude of the spin vector satisfies = S for 
the spin angular momentum S. The second term, which 
may be called the Zeeman term, is the scalar product 
of magnetization M = — and the magnetic 

field, — M ■ Bu,, where 5 = 2 and [is = efr/2mc. We 
note that the Weiss field (the internal magnetic field) is 
Bweiss — i^bJabS/ gfJ-B — 240 T, which is much larger 
than that of microwaves, where ns =4.1 is the effective 
magneton number of magnetite and Jab is of a few meV. 
We also note that the individual spin interaction energy 
JabS 2 = 0.016 eV is comparable to the thermal energy 
at room temperature 0.026 eV. Thus, thermal effects are 
significant in the sintering process for which temperatures 
are much above 300 K. 

To obtain a thermally equilibrated state, we minimize 
the internal energy of the spin system Eq.(l) at a given 
temperature using the Monte Carlo method with the 
Metropolis criterion. A trial random rotation is exerted 
on one of the spins in the n-th step: the trial is accepted 
if the internal energy decreases 6U = U n — U 71 ^ 1 < 0, 
or if the energy increase satisfies cxp(— SU/ksT) > e, 
where e is a random number uniformly generated in the 
(0,1) interval, /c^ is the Boltzmann constant and T is 
temperature; otherwise the trial is rejected. 

To proceed further, we have to specify the crystal 
structure and the exchange interaction coefficients. For 
magnetite, the unit cell is a cubic box with the sides 
of 0.8396 nm at room temperature^ and contains 24 
irons: 16 of them are Fe 3+ (3d 5 ,S = 5/2) and occupy 
all the tetrahedral sites (A-site, 8 positions) and half 
of the octahedral sites (B-site). The rest of irons are 
8 Fe 2+ (3d 6 ,S = 2), and are located at B-sites. The 
exchange interaction coefficients are all negative and sat- 
isfy | Jab | > \Jaa\,\Jbb\^- We assume Jab = -4.0 
meV and Jaa = Jbb = —0.3 meV, where the former 
is larger in magnitude than the theoretically estimated 
value 2.5 meV— to reproduce T c in our model. 
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FIG. 1: Temperature dependence of the magnetization along 
the c-axis (z-axis), M z — — '^2 i gfJ,BS Zt i, calculated by Monte 
Carlo simulation using the Heisenberg model for magnetite 
described in the text. The magnetization is normalized by 
Mo = NfiB , where N is the number of Fe ions. The magneti- 
zation vanishes above the Curie temperature T c of magnetite 
858K. 



III. RESULTS 

Using the procedures described in Sec. II, we first cal- 
culate the equilibrated state without an applied magnetic 
field for the periodic crystal of magnetite with 3x3x3 
unit cells. Below T c , spins are ordered along the c-axis 
due to the exchange interactions, with the spins at A-sites 
and those at B-sites oriented oppositely. The temper- 
ature dependence of magnetization calculated with our 
model is shown in Fig.l. The validity of our model is 
confirmed by the numerical result that the spontaneous 
magnetization decreases monotonically with temperature 
and vanishes above T c , as expected. Our calculated value 
M c /Mo = 1.1 at 300 K roughly agrees with the experi- 
mental value M s /M = 1.3±£, where M = N^ B with N 
the number of irons in the system. 

Next, we apply a slowly varying magnetic field B m to 
the equilibrium obtained above, and perform the Monte 
Carlo simulation. The varying speed of the magnetic field 
is chosen to be slow enough for spins to relax to the state 
of an instantaneous magnetic field. The amplitude of the 
magnetic field is 1 T unless otherwise specified, which 
is still by two orders of magnitude less than the Weiss 
field. Spins relax in orientation under the magnetic field 
and the internal energy U(B W ) depends on the magni- 
tude and orientation of B w . We evaluate the difference 
between the maximum and minimum of the internal en- 
ergy during a period of the magnetic field change, 
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Here, the superscript (mc) stands for "Monte Carlo" 
simulation. The energy AU that is to be released irre- 
versibly to lattice atoms by dissipation process, possibly 
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FIG. 2: Model calculation for the magnetic field effects on 
magnetite at 700 K. The changes in (a) the z-component of 
magnetization M z = — ^ g\iBS z ,i, (b) the internal energy U 
per Fe (solid line), and the contribution of the Zeeman term 
to the internal energy (dashed line, for which the zero baseline 
is shifted downward by 7!7o) are shown against (c) the applied 
magnetic field B z in the Monte Carlo simulation steps. The 
z-axis is parallel to the c-axis. The normalization factors are 
as follows: M = N/i B , U = 10~ 21 J, and B = 1 T. 



by spin-lattice interactions^, is obtained by subtracting 
from A[/( mc ) the reversible energy change, AU^ rev ^ = 
The latter is obtained by solving the 
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spin dynamics equation without dissipation, 
dsi/dt = y^(2Jij-/fi)Bj x Sj - (gus/tysi x 



(3) 



The energy difference AU can be used as the index of 
heating by the microwave magnetic field. This is verified 
later by agreement of the estimation based on the imag- 
inary part of the magnetic susceptibility with the AU 
here and also the experimental value of the heating rate. 

Figure 2 shows the changes in (a) the average magne- 
tization in the z direction (c-axis), (b) the total internal 
energy (solid line) and the contribution of the Zeeman 
term (dashed line) at 700 K against (c) the alternat- 
ing magnetic field which is parallel to the c-axis. In the 
present case, the magnetization M z stays negative along 
the z-axis, and the Zeeman term — M • B„, takes positive 
or negative value according to whether B„, is parallel or 
anti-parallel to the c-axis. In the phase where the applied 
magnetic field becomes more negative, the magnitude of 
magnetization increases due to alignment of spins paral- 
lel to the c-axis, and the internal energy becomes min- 
imal in this phase of the magnetic field. We note that 
the large change in the internal energy occurs through 
the exchange interactions because the energy associated 
with the Zeeman term is small and reversible, with the 
former by a factor of JabS/ gfJ-BB w (3> 1) larger than 
the latter, as depicted in Fig. 2(b). In fact, the contribu- 
tion of the Zeeman term is small and roughly the same 
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FIG. 3: The distribution functions of spins along the x di- 
rection (leftward panels) and z direction (rightward panels), 
G is the angle of the spin either with the x- or z-axis. The 
distribution functions in the upper panels correspond to Fe 3+ 
(the shaded area denotes the states of spins at A-sites) and 
those in the lower panels correspond to Fe 2+ . 
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FIG. 4: Temperature dependence of the internal energy 
change (per Fe) during a period of an applied magnetic field of 
B w = 1 T (U = 10" 21 J). The filled circles and squares cor- 
respond, respectively, to the irreversible energy difference A U 
and the change in the Zeeman term A(M • B w ) for magnetite 
with the magnetic field parallel to the c-axis; crosses corre- 
spond to A(7^ rel '' obtained by the spin dynamics of Eq.(3) 
which is reversible. The open circles and squares correspond 
to AU and the change in the Zeeman term, respectively, when 
the magnetic field is applied parallel to the a-axis. The verti- 
cal line denotes T c of magnetite. 



as that obtained with the dissipationless spin dynamics 
Eq.(3). 

The distribution functions at the 1.25 x 10 6 th Monte 
Carlo step in Fig. 2 are shown in Fig. 3. The x and z di- 
rections are taken along the a- and c-axes, respectively. 
In the x direction, the spins form centered Boltzmann 
distributions and no magnetization occurs along this di- 
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FIG. 5: Temperature dependence of the internal energy 
change (per Fe) AU for hematite is shown by triangles, and 
that of magnetite with the magnetic field parallel to the c- 
axis is shown by filled circles as reference, for B w = 1 T 
(Uo = 10 -21 J). The imaginary part of the magnetic sus- 
ceptibility obtained by dissipative spin dynamics Eq.(4) for 
to = 1 ns is shown by double circles (in mol -1 ). The vertical 
line denotes T c of magnetite. 



rection. (The shaded areas correspond to the Fe at 
A-sites.) Fcrrimagnctization occurs along the z direc- 
tion. The Fe 3+ spins at A-sites are anti-parallel to the 
z-axis and those of Fe 2+ and Fe 3+ at B-sites are paral- 
lel to it. The change in the orientations in response to 
the microwave magnetic field is subtle, yet it gives rise to 
a maximum in the magnitude of magnetization periodi- 
cally, as shown in Fig. 2 (a). 

Figure 4 shows the calculated temperature dependence 
of the internal energy difference AU = A[/( roc ) - AU^ 
when the microwave magnetic field is parallel either to 
the c-axis (filled circles) or to the a-axis of magnetite 
(open circles). The difference in the energy becomes 
largest when the polarization of the magnetic field is 
parallel to the c-axis, and it increases linearly with tem- 
perature up to T c . The change in the reversible energy 
A.U( rev ) anc i that in the Zeeman term for the former case 
are shown by crosses and filled squares, respectively The 
energy change AU^ rev ^ is similar to that in the sponta- 
neous magnetization which decreases with temperature 
and vanishes above T c . The change in the Zeeman term 
is almost the same as AU^ rev ^ but is finite for T > T c . 
The Zeeman term in the Monte Carlo simulation re- 
mains finite because an induced magnetization appears 
synchronously with and along the magnetic field inde- 
pendently of its polarization in the paramagnetic regime. 
The major contribution to the internal energy difference 
AU is attributed to the exchange interactions since the 
change in the Zeeman term is reversible and thus almost 
subtracted in AU. The temperature dependences shown 
in Fig. 4 arc in excellent agreement with the sintering ex- 
periments in the microwave magnetic field 3 -^, where the 
heating of magnetite was enhanced at 300 — 600° C and 



FIG. 6: The relaxation time td dependence of the imaginary 
part of magnetic susceptibility \ ( m mol -1 ) for the 2.5GHz 
microwave and temperature 700K. 
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FIG. 7: The microwave frequency uj dependence of the imag- 
inary part of magnetic susceptibility \ ( m mol -1 ) for the re- 
laxation time td = 1 ns and temperature 700K. 

continued to much above T c . 

An argument that involves time scales as well as the 
temperature dependence of heating is made possible by 
solving the dissipative spin dynamics. Here, the term 
— (si — Sio)/rD is added to the righthand side of Eq.(3). 
Namely, we solve 

dsi/dt = y~](2Jjj /fj)8j x sj - (gfi B /fi)s t x B w 

-(Si -Sio)/r D . (4) 

A set of Sio constitutes the equilibrium spin distribution 
function for a given temperature and magnetic field un- 
der the Monte Carlo calculation. The imaginary part of 
the magnetic susceptibility Xi is obtained from the linear 
response of magnetization against the applied alternat- 
ing magnetic field. Xi gives the heating rate dT/dt cx Xij 
and is shown for magnetite with double circles in Fig. 5. It 
peaks around 800 K and is consistent with AU obtained 
with the energy principle of the Monte Carlo calculation. 
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FIG. 8: Magnetic field dependence of the internal energy 
change in magnetite per wave period for the 2.5GHz mi- 
crowave and temperature 700 K. The dotted line corresponds 
to B 2 . 

The dependence of the imaginary part of magnetic sus- 
ceptibility Xi on the relaxation time td is shown in Fig. 6 
for the 2.5 GHz microwave. We see that Xi 1S inversely 
proportional to td in the r = 2it /ld < td range. Also, the 
dependence of Xi on the microwave frequency is shown 
in Fig. 7 for td = 1 ns. The imaginary part of mag- 
netic susceptibility peaks around 2 GHz. This agrees 
with the experimentally deduced magnetic permeability 
for magnetite^. In our calculation, the relaxation time 
is assumed to be constant irrespectively of temperature 
or microwave frequency. However, if the relaxation time 
becomes small with temperature rise, the heating rate in- 
creases as r^ 1 as shown in Fig. 5, resulting in more rapid 
heating at elevated temperatures. 

From these data, one obtains the heating rate dT/dt ~ 
(l/2)u>XiB^ nw /c p = 300 K/s for the case of td = 1 ns 
and the microwave field B mw = 150 G, where the heat 
capacity c p is 210 J/K mol at 600K (the wave period is 
400 ps for the 2.45GHz microwave). This is large enough 
to account for the experimental value (dT/dt) ex = 250 
K/s in the microwave sintering of magnetite^. 

The dependence of the internal energy change in mag- 
netite per wave period on the applied magnetic field is 
shown in Fig. 8. The magnetic field amplitude ranges 
from 0.125 T to 1.0 T. It is well fitted by the quadratic 
law B 2 , therefore, AU thus the heating rate is propor- 
tional to the microwave power. 

Hematite has a different crystal structure from 
magnetite^ and has weak spontaneous magnetization. 
The calculated change in the internal energy AU is plot- 
ted for various temperatures by triangles in Fig. 5. By 
comparison of AU with the case of magnetite, the re- 
sponse of hematite to the alternating magnetic field is 
considerably weak. This is consistent with the exper- 
imental fact that only the domain of magnetite in a 
magnetite-hematite composite powder was heated by the 



magnetic field of microwaves^. 

The observation that heating of titanium oxide occurs 
only when oxygen defects are present^, Ti02-:c (x > 0), 
is explained in a, similar fashion. Titanium in titanium 
oxide Ti02 is a Ti 4+ ion with 3d electron configura- 
tion, which has no electrons in the 3d shell. Thus, it has 
no spontaneous magnetization and should not respond 
to microwave magnetic field, similarly with hematite. 
However, when there are oxygen defects, the trivalent 
titanium ions Ti 3+ appear which have 3d 1 electron con- 
figuration. The spins of these 3d electrons respond to 
the magnetic field of microwaves, and absorb microwave 
energy which leads to observed heating under the mi- 
crowave magnetic field. 

IV. SUMMARY 

In this paper, we have showed theoretically the mech- 
anism of the rapid and selective sintering of magnetic 
metal oxide particles by the magnetic field of microwaves. 
We adopt the Heisenberg model, and perform both 
Monte Carlo calculation and dissipative spin dynamics 
simulations. The heating occurs due to the response of 
magnetization to microwaves, which originates from elec- 
tron spins residing in the unfilled 3d shell. Their non- 
resonant response causes a large change in the internal 
energy through the exchange interaction between spins. 
It persists above the Curie temperature T c because each 
electron spin is able to respond to the alternating mag- 
netic field of microwaves even above T c . This energy 
change will then be dissipated to lattices and contribute 
to heating. 

Hematite FC2O3 which has only weak spontaneous 
magnetization shows much less response to microwaves 
than magnetite. Also, the heating of titanium oxide hav- 
ing oxygen defects T1O2-X (x^0) by the microwave mag- 
netic field is explained by our theory in terms of intrinsic 
(spontaneous) magnetization. 

The imaginary part of the magnetic susceptibility Im% 
obtained by solving dissipative spin dynamics agrees well 
with the heating results by the Monte Carlo calculation. 
We have also presented the dependences of the heating 
rate on the frequency of microwaves and on the spin re- 
laxation time. These results well account for the large 
heating rate of magnetic metal oxide by the microwave 
magnetic field in the sintering experiments. 
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